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Glossary 
Cluster 

Clusters are sets of neighbouring sites of the same type. 
Ising model 

Each site carries a magnetic dipole which points up or down; neighbouring 
dipoles "want" to be parallel. 

Opinion dynamics 

How do people change opinions? Simulations usually ignore all details of 
the brain and represent the opinion by one or several numbers which can be 
changed due to contact with others. 

Schelling model 

People belonging do different groups may produce segretated neighbour- 
hood just by their personal preferences, not by outside force. 

Sociophysics 

Application of methods from (mostly statistical) physics to human rela- 
tions; can be traced centuries backwards. 
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Figure 1: Ising model after 20 Glauber kinetic steps per site on a 500 x 500 
square lattice at fc^T/J = 2. We start from a random distribution of equally 
many black and white sites. Figs.l, 2, 3 omitted for arXiv; their small are 
not readable on some computers and show the standard phase separation or 
clustering of Ising model.. 



Figure 2: As Fig.l but after 2000 instead of 20 iterations. 

1 Definition and Introduction 

The application of concepts from the natural sciences to social sciences, partly 
to be reviewed here, is at least 25 centuries old. Then the Greek philosopher 
Empedokles stated (according to J. Mimkes) than humans are like liquids: 
Some mix easily like wine and water, and others like oil and water refuse to 
mix. We start with the Schelling model of 1971 which implemented this idea, 
and its criticism. Then we will review opinion dynamics in large populations, 
summarising only shortly other aspects like self-organisation of hierarchies 
or competition between human languages. 

Humans do not like to be treated like a number, and indeed the human 
brain is much more complex than a binary variable (called "spin" by physi- 
cists) which is either +1 or -1. We do not deal here with the psychological 
processes of an individual but with mass psychology, and this author learned 
half a century ago in school that mass psychology is different from individual 
psychology: The law of large numbers averages out over individual fluctua- 
tions and makes general trends more clearly visible. Thus what we call today 
statistical physics plays a useful rule, and social scientists [U [2] have applied 
it, without knowing then that they dealt with an Ising model of ferro magnets. 

The astronomer Halley, best known through his comet, tried to establish 
mortality tables already three centuries ago. Of course, the time of death 
of one given individual is usually difficult to predict but averaged over mil- 
lions of people the statistical offices of many countries prepare regularly life 
tables which tell us how probable it is for a newborn child to live up to x 
years, provided there are no changes of the mortalities in the coming decades. 
Insurance for automobiles is another example: We do not want to produce 
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accidents, but we know that they happen, and take precautions for their fi- 
nancial consequences. Thus the whole insurance industry is based on treating 
humans like numbers, ignoring their individuality. 

Finally, human opinions are often fluctuating and ill-defined, but never- 
theless in elections people cast one choice, out of a limited number of choices. 
And election results belong to those social data for which we have lots of ac- 
curate numbers, based on large populations. 

Thus it is not at all the merit (or ignorance) of physicists which treats 
humans like numbers; this method has a very long tradition and is an indis- 
pensable part of modern life. 

2 Schelling Model 
2.1 Ising simulations 

Following (but not citing) Empedokles, the later economics Nobel laureate 
Schelling [1] asked whether the racial segregation in American cities can 
emerge from intrinsic behaviour of the individual people, instead of or in 
addition to extrinsic reasons like discrimination, rent differences, etc. In 
particular, can "black" ghettos in the predominantly "white" USA arise just 
because people prefer to have neighbours of their own group over neighbours 
from the other group? In many other countries we find many other types 
of residential segregation, based on religion, ethnicity, .... In physics, such 
a process is easily simulated through the two-dimensional Ising model, as 
shown in Figs.l and 2. 

In this Ising model, each site on a square lattice carries a variable Si = ±1, 
and each pair < i,k > of nearest neighbours produces an "energy" —JSiSk 
with some proportionality constant J. The total energy E (= total unhap- 
piness) is the sum of this pair energy over all neighbour pairs of the lattice. 
In physics, different distributions of the "spins" Si are realised with a prob- 
abihty proportional to exp^—E/ksT) where T is the absolute temperature 
and ks the Boltzmann constant. There is no need to worry about values for 
T, ks, J since the only relevant quantity is the ratio ksT/J, taken as 2 in 
these pictures. The "Glauber" kinetics is simulated on the computer by flip- 
ping a spin if and only if a random number between and 1 is smaller than 
the probability exp{—AE/kBT)/[l + exp{—AE / ksT]. The Fortran program 
contains less than 40 lines and takes a few seconds: 
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parameter (L=500 , Lmax= (L+2) *L) 
dimension is (Lmax) , iex(9) 
byte is 

data t, max, ibm/2. 00, 2000,1/ 

print *, L,max,ibm,t 

Lpl=L+l 

L2pL=L*L+L 

do 1 i=l,Lmax 

is(i)=-l 

ibm=ibm* 16807 

1 if(ibm.gt.O) is(i)=l 
do 2 ie=l,9 

ex=exp(-2* (ie-5) /t) 

2 iex(ie)=(2.0*ex/(1.0+ex) - 1 . 0) *2147483647 
ibm=2*ibm+l 

do 3 mc=l,max 
do 4 i=Lpl,L2pL 

ie=5+is(i)*(is(i-l)+is(i+l)+is(i-L)+is(i+L)) 
ibm=ibm* 16807 

4 if (ibm. It . iex(ie) ) is(i)=-is(i) 
mag=0 

do 6 i=Lpl,L2pL 
6 mag=mag+is (i) 

3 if (mc.eq. (mc/100000)*100000) print *, mc.mag 
do 5 i=Lpl,L2pL 

if (is(i) .ne . 1) goto 5 
iy = i/L 
ix=i-L*iy 
print *, ix, iy 

5 continue 
stop 
end 

Such models and programs are taught in courses on computational or 
theoretical physics all over the world; the model was published in 1925. If 
in the above flipping probability the denominator is omitted one gets the 
Metropolis kinetics. If instead of flipping one spin we exchange two opposite 
spins, we get the Kawasaki dynamics. For Glauber or Metropohs, after very 
long times (measured by the number of sweeps through the lattice) one of 
the two possibilities dominates at the end, if T is not larger that the critical 
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temperature T^, with 2J/kBTf. = ln(l + \/2) ~ 0.44 known since 1940. For 
Kawasaki dynamics the fraction of black sites remains constant, and we get 
two large domains. For higher temperatures above Tc only small clusters and 
no large domains are formed, Fig.3. 



Figure 3: As Fig. 2 but at ksT/J = 3 instead of 2. Only small clusters and 
no large domains are formed. After 200 and 20,000 iterations the pictures 
look similar to this one made after 2000 iterations. 

In this Ising model, two neighbouring spins have due to their interaction 
—JSiSk a higher probability to belong to the same group than to belong to 
the two different groups. If the difference between these two probabilities is 
large enough, T < T^, domain sizes can grow to infinity in an infinite lattice. 
Figs. 1 and 2, while only small clusters are formed in Fig.3 for smaller 
differences in the probabilities, T > Tc. That these probabilities, controlled 
through —J/ksT, lead to these different regimes, separated by a sharp phase 
transition at T = Tc, is not obvious from the definition of the interaction 
JSiSk, took physicists many years to find, and is typical of complex systems. 

The social meaning of temperature T is not what we hear in the weather 
reports but an overall approximation for all the more or less random events 
which influence our decisions but are not explicitely included in the model. 
For residential segregation the model only counts how many neighbours of 
which group one has. But not all people of one group are alike, housing in 
different parts of a city costs different amounts of money, some parts are 
more beautiful then others, and job hunting may force us into a temporary 
residence of a new city which does not conform to our wishes. In this way, a 
positive temperature allows for rare moves which increase the energy, i.e. we 
move to a new residence where the neighbourhood composition along makes 
us less happy. At zero temperature, the Ising model does not properly order 
into one or two "infinite" domains. 

2.2 Schelling's version and later improvements 

Schelling pj avoided probabilistic rules and thus counted neighbours = ±1 
at zero temperature. Then it does not matter if all neighbours or only a 
majority of them belong to the own group. Thus people are defined as happy 
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if at least half of the neighbours belong to the own group, and as unhappy 
otherwise (i.e. if the majority belong to the opposite group). Unhappy 
people move to the nearest place where they are happy. Since Schelling 
moved only one person (or family) at a time, and made no exchange of two 
people simultaneously as in Kawasaki kinetics, he introduces a large fraction 
of empty residences. Thus at each step, one unhappy person or family moves 
into the closest vacancy where hfe would be happy. 

This model, and also many variants [H [3], fail to give large domains; 
only small clusters are seen. In reality, Harlem in Manhattan (New York), 
is not a cluster of a few houses but extends over many square kilometers. 
Thus the original version does not give the desired results. Large domain 
are formed if people also change residences if this brings no improvement 
|3] (hardly a realistic assumption) or at a finite temperature |5]. The latter 
paper also gives some alternatives to the Schelling model which also allow 
for large domains, and a simple example of s finite cluster where everybody 
is "happy" and which therefore never grows or dissolves on its own "will" . 

Much earlier and simpler is the zero-temperature version of Jones [2] who 
at each iteration removes a random fraction of the people and fills the vacan- 
cies with people who are there happy in the Schelling sense. This randomness, 
just like the finite temperature, leads to large domains as desired. Neither 
physicists nor social scientists have taken much note of [2]. The history of 
the Schelling model is an example how the lack of communication between 
disciplines has hampered progress in research, even very recently [H [S] . Only 
computational statistical physicists know everything. (Jones in footnote 4 of 
p] also mentions a probabilistic version closer to the Ising model.) 

For finite temperatures, [_5j follows the above Glauber dynamics, but in- 
stead of an energy E uses a variable which is or 1 depending on the hap- 
piness of the residents. Moving from one place to the other then depends 
exponentially on the ratio of this variable to ksT, instead of on the ratio 
AE/ksT. Many variants are possible, e.g. in the treatment of neutral cases 
[H [6] where the number of neighbours of both groups is exactly the same. 

But we are on safer grounds and can use decades of physics research if 
we use the normal Ising model, or its generalisation to Q different groups, 
the Q-state Potts model. Then [71 18] implemented a suggestion of Weidlich 
[9j that people slowly learn to live together with neighbours from the other 
group. Thus T not only takes into account the various accidents from outside 
the model, but also measures the tolerance: The higher T is the more are 
people willing to live in neighbourhoods of the other group. In the limit 
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T = oo the neighbours would not matter at all, for intermediate T, Fig. 3 
showed small clusters but no large domains, and for low T the domains grow 
to infinite sizes on an infinite lattice. The learning suggested by Weidlich 
thus means that this parameter T (= temperature or tolerance) no longer is 
kept constant but slowly increases. 

For an Ising model, [7| showed how an initial large domain dissolves if the 
temperature is slowly increased from below to above Tg. More realistically, 
for five (instead of only two) different groups in a modified 5-state Potts 
model, [8] increased T from low to high values and showed that with a slow 
increase one has appreciable domain formation during intermediate times, 
while with a fast increase this segregation is mostly avoided. Fig. 4. 
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Figure 4: Amount of neighbours of the same type in a Potts model of five 
groups, normalised to unity for the initial random distribution. The temper- 
ature or tolerance increases from low to high values, slowly in the top curves, 
and fast in the lower curves; the latter mostly avoid the segregation into dif- 
ferent group. (The lowest line holds for a constant high temperature.). From 

rs]. 
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3 Opinion Dynamics 



The following section describes several rules for simulated people to change 
opinions; each of these rules is applied again and again to these agents until 
some stationary or static state has been achieved. 

3.1 Ising model 

Also for human opinions, one could use the Ising model of the previous sec- 
tion [inilTT]; see also [9]. People can vote for or against the government or a 
new constitution, for one of two presidential candidates, or (using generalised 
models) for one out of Q > 2 different parties. Their neighbours on a lat- 
tice influence them in their vote, and in addition mass media may influence 
everybody in one direction. The latter effect can be modelled through an 
external "magnetic" field, eq.(lb) in "Phase transitions ..." by this author 
in this encyclopedia. No motion of people needs to be taken into account, 
and the complications of Kawasaki kinetics (exchange of two people with op- 
posing opinions) are not needed. Thus the Glauber program of the previous 
section still can be used, and we only refer here the old generalisation into 
the social impact model [T2| [13] and to a recent financial application |14j . 

3.2 Voter model 

Also quite old is the voter model [15]: Each person chooses between two 
opinions, by taking over the one of a randomly selected neighbour. One 
may rewrite this rule as stating that each person selects the opinion of the 
neighbourhood, with a probability proportional to the number n of neigh- 
bours having that opinion. Thus in contrast to the Ising model where the 
probabilities depend exponentially on ra, now they depend linearly on n. A 
final equilibrium (absorbing fixed point) is reached if everybody shares the 
same opinion. The deviations from that final state can be measured by the 
magnetisation (difference between the numbers for the two opinions) or en- 
ergy (average number of neighbours having the opposite opinion). The time 
needed to reach the consensus increases with a power of the lattice size, and 
the exponent depends on the dimensionality. A nice and short review of the 
voter model, also on various networks, is given by the Majorca group |16j . 
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3.3 Axelrod model 



Axelrod [T7] wondered how different opinions or cultures may coexists even if 
people tend to become more alike in their beliefs. Looking at the above Ising 
Figs. 1, 2, 3, we see that due to finite time and/or finite temperature such 
coexistence of two opposing opinions is possible. But Axelrod generalised it 
not only to Q > 2 different possible opinions as in the Potts model of the 
"Schelling" section, but also to F different questions. People may have one 
set of opinions on which political party they want to vote for, another set 
about what is the best football team, a third about recent cinema films, etc. 
This allows for different opinion sets on all F questions ("features"). Of 
course, one could could generalise this model to the case where the number 
Qf of possible choices is different for the different features /, allowing then 
for Uj^iQf instead of simply different sets of opinions. 

Another aspect of the model takes into account that people prefer to 
talk to, or to make political coalitions with, others with whom they share 
many opinions. Thus the probability of one person to take over the different 
opinion of a neighbour is proportional to the number of features on which 
their opinions already agree. In the next subsection we will use a similar 
concept under the name of bounded confidence. 

Whether a total consensus ( "globalisation" ) is reached or multiculturality 
persists depends on parameters: Small Q lead to consensus. Again, the 
Majorca group [TB] reviewed the many follow-up papers on this Axelrod 
model. 

3.4 Sznajd, Krause-Hegselmann and DefFuant models 

Much of the opinion dynamics research since 2000 centred on three different 
models S, KH and D, originally invented independently around that year: 
Sznajd [H] (S), Krause-Hegselmann [19] (KH) and Deffuant et al |20] (D). 
They were also called missionaries, opportunists and negotiators by some 
computational physicists [39j. 

The S model is closest to the earlier models since it allows for Q discrete 
opinions, while KH and D use real opinions, e.g. between zero and one. S 
happens on a lattice or network while for KH and D everybody may interact 
with everybody. In the most widespread S version a pair of neighbouring 
sites on a square lattice convinces its six neighbours of its opinion, if and only 
if the two opinions of the pair agree [21]; governments and parties usually 
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lose support if their internal opinion differences make it to the headlines. 
For KH, the new opinion of a person is the arithmetic average over the 
opinion of the whole population. For D, each person selects randomly another 
person and then both move in their opinion towards each other by an amount 
proportional to their opinion difference. 

In all three cases, "bounded confidence" applies: The KH agents average 
only over those people who differ from their own opinion by less than e, and 
the D agents only select negotiation partners differing by less than e from 
their own opinions. In both models < e < 1 is a fixed parameter. For S 
agents with Q = 2 such a rule makes no sense, but for Q > 2 one can modify 
the convincing rule such that only neighbours differing by at most ±1 from 
the pair opinion adopt the pair opinion. Thus 1/Q for S plays the role of 
e for D and KH. A rule similar to this bounded confidence was mentioned 
above for the Axelrod model [17j . 

Inspite of the differences in their definitions, the results are quite similar 
for S, KH and D. For large e or Q < 3 a complete consensus is usually 
reached; for small e or Q > A different opinions may coexist forever. In 
addition to computer simulations, also analytical calculations were made 
[22l [23] which agree with many aspects of the simulations. More results, 
also for opinions on more than one feature and agents sitting on scale-free 
networks, are summarised by us in [39]. 




One particular application is shown in Fig. 5: Various Brazilian election 
results for candidates in city councils showed great similarity if the number of 
candidates getting a given number v of votes is plotted against this v. Putting 
the S model with Q = 1000 candidates onto a scale-free network instead of 
a square lattice, excellent agreement of simulation and reality was found 
after the numbers were scaled by suitable factors. Also Indian elections were 
simulated this way [2S], while some exceptions may also exist (S. Fortunato, 
unpublished). It would be nice to apply other opinion dynamics models to 
the same election problem. As usual in statistical physics, these studies can 
predict and simulate the shape of the distributions but not the winner in 
a specific election, just as we can predict the pressure of the air molecules 
around us but not where which molecule will be one minute from now. 

3.5 Galam conservatism 

Galam has published since many years theoretical models which may explain 
why reforms are very difficult and why a minority can stay in power. Usu- 
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Figure 5: Brazilian elections (x) and simulations of Sznajd model on Barabasi 
- Albert networks (+); from |24j . 

ally these models are solvable analytically and assume that the population is 
divided into small groups of people which to the outside are represented by 
one person who follows the majority wish of the group. Several of these rep- 
resentatives form a supergroup, and this supergroup again decides according 
to the majority of the representatives in it. In this way an "infinite" hierar- 
chy of people, groups, supergroups etc can be built. In the case of equally 
many voting for one choice as for the opposite choice, within one unit, that 
unit votes for the status quo. Starting with everybody having opinion -1, a 
very large majority of people must switch to opinion +1 before the top of 
the hierarchy finally also changes opinion |26j . We refer to our book ^39j for 
a summary of more recent Galam papers. 
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4 Languages, Hierarchies and Football 



4.1 Language competition 

Darwinian survival of the fittest is established biology, but similar concepts 
can be applied to human languages, bridging the gap to opinion dynamics. 
There are now thousands of different languages, and their "size" is the num- 
ber of native speakers of that language. The size distribution extends from 
1 (on the verge of extinction) to 10^ (Mandarin Chinese). The grammar of 
a language |27j can be characterised by F features each of which can have 
Q different values, just as in Axelrod's model explained above. Features 
can change spontaneously or be taken over from a (neighbouring) language; 
speakers of a small language give it up and learn a widespread language (as 
done with physics research publications since 1945); people migrate to other 
places and bring their language with them. All these processes can lead to 
the extinction of existing languages and the creation of new ones (by the 
branching of one language into several daughter languages.) The present 
language size distribution is roughly log-normal, with an enhancement at 
small sizes [2S] • Similar languages form families, and the size distribution of 
families is a power law at intermediate sizes [22] (where the size is now the 
number of different languages) belonging to that family). 

Various computer simulations of this language competition have been 
made, mostly since 2003 and reviewed recently [30]; see also [311 [32]. We only 
mention Fig. 6 from a modified Viviane model which agrees well with the real 
language size distribution. For language families, the empirical statistics is 
worse [29] and the models work less well |30j . 

4.2 Self-organisation of social hierarchies 

The elites of all countries and all times always had excellent reasons why they 
should be on top and others on the bottom. This holds even when the United 
Nations criticise the school system as violating human rights. In contrast, 
the Bonabeau model [M] explains social hierarchies as purely accidental, 
without any merit. People are put on a lattice, occupying a fraction p of all 
lattice sites and having an initial score of zero. Then they move randomly 
to neighbouring sites, and whenever one person wants to move into the site 
occupied by another person, a fight erupts. The winner takes the contested 
site, the loser moves into (or stays at) the other site. Also, the winner adds 
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Figure 6: Simulated language size distribution on a 20, 000 x 20, 000 square 
lattice using a modified Viviane model [33] . 

one point and the loser subtracts one point in its score, and in the future the 
agents with a positive score have a higher probability to win, those with a 
negative score have a lower probability to win. Slowly the history is forgotten, 
by reducing the score at each time step by, say, ten percent. 

With some suitable feedback between the distribution of scores and the 
probability to win, a phase transition was simulated such that for p above 
some critical concentration, the standard deviation in the scores becomes 
positive for long times and large populations. For p < pc it fluctuates near 
zero, that means everybody remains near a 50 percent chance to win. So, 
just by accident at a high population density some people rise to the top, 
and others fall to the bottom. However, the people on top (bottom) are 
not always the same; only the differences between top and bottom, not the 
people, remain the same. [551 ESI EZl are some of the more recent references 
in this field. 
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4.3 Football 

Football (= soccer) is the world's most popular spectator sport, though in 
the author's city it is more a frustration. Randomness surely plays a role 
and makes it attractive. Can we explain all results just by chance, in the 
spirit of Bonabeau hierarchies? Assuming a constant probability to make a 
goal within one minute, the distribution of goals and victories is more narrow 
than in reality. If instead we assume that this probability varies from team 
to team, still no good agreement is found. Good agreement with reality is 
obtained only if correlations are taken into account [38], in the sense that a 
goal makes the scoring team happy, shocks the opposing team, and thus with 
an enhanced probability leads to another goal for the scoring team. Thus if 
we lose it is not just bad luck; it is also the referee's fault. 

5 Future Directions 

The Schelling model of section 2 is not the only case of missed opportuni- 
ties because of a lack of cooperation between social sciences on the one side 
and physics, mathematics or computer science on the other side. The two 
books cited at the end |39j were written without the authors of one book 
knowing of the preparation of the other book. One group of authors works 
in physics departments; none of the other group lists physics as institutional 
address. Nevertheless the two books show strong overlap in fields and meth- 
ods covered, but little overlap in the literature cited. More interdisciplinary 
cooperation would help. 
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